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We obtain a magnetically charged regular black hole in general relativity. The source to 
the Einstein field equations is nonlinear electrodynamic field in a physically reasonable model 
of nonlinear electrodynamics (NED). “Physically” here means the NED model is constructed 
on the basis of three conditions: the Maxwell asymptotic in the weak electromagnetic held 
limit; the presence of vacuum birefringence phenomenon; and satisfying the weak energy 
condition (WEC). In addition, we analyze the thermodynamic properties of the regular black 
hole in two ways. According to the usual black hole thermodynamics, we calculate the heat 
capacity at constant charge, from which we know the smaller black hole is more stable. We 
also employ the horizon thermodynamics to discuss the thermodynamic quantities, especially 
the heat capacity at constant pressure. 

PACS numbers: 04.70.Dy 


I. INTRODUCTION 

The first well-known model of nonlinear electrodynamics(NED) is the Born-Infeld theory(BI), 
which is proposed to obtain a finite electron self-energy[l]. Heisenberg and Euler found that 
due to the presence of virtual charged particles the one-loop quantum correction in quantum 
electrodynamics will give nonlinear contribution[2]. Not only that, the virtual particles will result 
in “ polarization of the vacuum ”. In this case, the vacuum behaves like a polarizable continuum and 
should exhibit the phenomenon of birefringence[3-5]. These effects can be observed in experiments 
such as PVLAS[6] and BMV[7] and the experimental results can put some restrictions on the 
parameters introduced in the NED models. Vacuum birefringence is a nonlinear effect. Therefore, 
we expect that in a physically reasonable model of NED, there should exist the effect of vacuum 
birefringence. In this sense, the usual BI theory is not a physically allowable NED model because 
of the absence of vacuum birefringence. However, in the generalized BI model with two parameters 
the vacuum birefringence may exist [8]. The nonlinear effects in electrodynamics are significant only 
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for strong electromagnetic fields. In weak field case, NED should return back to the conventional 
Maxwell theory, just like the weak field approximation of general relativity gives the Newtonian 
mechanics. Thus, we expect that this is another requirement for a physically reliable NED model. 

Recently, due to the natural emergence in string theories [9-11], NED is coming into view in 
gravitational theories as source. It is well-known that by minimal coupling of Maxwell electromag¬ 
netic fields to gravity, the Reissner-Nordstrom(RN) black hole can be obtained. Similarly, more 
interesting black hole solutions can be derived through minimal coupling to gravity of nonlinear 
electromagnetic fields, such as BI and Bl-like electromagnetic fields[12, 13], logarithmic electro¬ 
magnetic field[14], power electromagnetic field[15] and exponential electromagnetic field[16]. 

In this letter, we shall investigate regular black hole solution derived in general relativity cou¬ 
pled to NED. The first example of a regular black hole was constructed by Bardeen in 1968[17]. 
Nearly thirty years later, Ayon-Beato, et al reobtained the Bardeen black hole by describing it as 
the gravitational field of a kind of nonlinear magnetic monopole[18]. Similarly, many other regular 
black holes can also be constructed by introducing nonlinear electromagnetic sources[19-23]. There 
are also other types of regular/nonsingular black holes with different origins. Eor more detailed de¬ 
scription of regular black holes, one can refer to the paper by Lemos et.al[24] and references therein. 
We should stress that there is a theorem which asserts that the existence of electrically charged, 
static, spherically symmetric solutions with a regular center is forbidden, while the existence of the 
solutions with magnetic charges is feasible, if the NED model contains the Maxwell theory as its 
weak approximations[22]. Thus, we only concern with the magnetically charged regular black hole 
with the special emphasis put on its thermodynamics and stability. 

The paper is arranged as follows: in the next section we simply introduce the NED model and its 
origin. In section 3 we will solve the Einstein field equations to obtain the regular black hole solution 
and analyze the geometric structure. In section 4 we verify that the nonlinear electromagnetic field 
in this NED model satisfy the WEC. Then we calculate the thermodynamic quantities, such as 
temperature, the heat capacity at constant charge and the heat capacity at constant pressure from 
which we can discuss the local stability of the regular black holes in section 5. We will make some 
concluding remarks in section 6. 
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II. THE NED MODEL 

Firstly we simply introduce the model of NED proposed by Kruglov[25], which can produce the 
vacuum birefringence phenomenon. The Lagrangian is given by 

where a is a dimensionless parameter and /3, 7 are parameters with the dimensions of [L^]. T = 
, Q = Ffj_^Ff^'' are two Lorentz invariants with and — 

is the electromagnetic field strength. 

It is shown that in a constant and uniform external magnetic field Bq, the indexes of refraction 
with different polarizations of electromagnetic waves are 


„ _ / + 1 )^ 


n± = 1 . 


This means that the electromagnetic waves with different polarizations have different velocities 
and thus the vacuum birefringence is present. However, one can easily see that this NED model 
cannot return back to the Maxwell theory in the weak approximation except for a = 0, because 

C ^ {1 + a)F + 0[F% (3) 

Therefore, below we shall consider a deformed NED model. 


2I3F+1 


7^2 

-F ■ 


Obviously, this model contains Maxwell theory as weak field approximation. Moreover, this model 
can still produce the vacuum birefringence phenomenon^. 


III. THE MAGNETICALLY CHARGED REGULAR BLACK HOLE 

We consider the NED in the framework of general relativity, with the action 

^ = -^1 [R - c{F, g)] , (5) 

where C{F,g) takes the form of Eq.(4). 

We only consider the static, spherically symmetric spacetime, in Schwarzschild gauge which can 
be written in the form 

ds‘^ = —f{r)dt^ + f{r)~^dr‘^ + r^dVL^. (6) 

^ Kruglov points out this fact for us. 
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The dynamical equation for the electromagnetic fields is 

= 0 , 

where Lp = dC/dF, Lq = dCjdQ, and the corresponding Bianchi identity is 

= 0 . 


( 7 ) 

( 8 ) 


According to the no-go theorem proposed by Bronnikov, to find the regular black hole solution, 
we will only consider the pure magnetic case[22]. In this case the ^-square term in Eq.(4) is zero. 
Thus the Lagrangian of electromagnetic fields now becomes 


C{F) 


F 

2^F + 1 ■ 


(9) 


In spherically symmetric case, involves a radial magnetic field T 23 and satisfies 


T23 = g'msin 6 ». 


( 10 ) 


where qm is the magnetic charge. For clarity and simplicity, we write Qm as q and take it to be 
positive below. Thus, F = 2F2sF^^ = 2g^/r^. 

From Fq.(9), one can easily seen that C{F) F, Lp ^ 1 when T" —)■ 0. This is the Maxwell 
asymptotic. Particularly, C{F) ^/‘2/3 = const when F —)• 00 , which means that near r = 0 the 
energy-momentum tensor T^i, is a constant. According to the discussion in [22], the spacetime near 
r = 0 is the de Sitter one and must be regular there. 

The metric function can be written in the form 


f{r) = 1 - 


2m{r) 


Substituting it into the Finstein field equations, one can obtain 


^0 ^1 2 m'(r) ^0 

Lr 0 — Of 1 — 2 — 0 — 1 — 


Integrating the above equation, we have 

1 

m{r) = M — - r'^Cdr, 
^ Jr 

where M is an integration constant and is chosen to satisfy m( 
the field equations of ( 2 ) and ( 3 ) components automatically. 
The mass function takes the form 

[in ~ 


-C/2. 


( 11 ) 


( 12 ) 


00 = 


(13) 

M. The solution (13) satisfies 


m(r) = 


- in-^- 

16^ [ 2^q + 2^^r + r2 


-|- 2 tctn ^ I 1 “h 




— 2 tan ( 1 — 




(14) 



5 


To obtain a black hole solution which is regular at r = 0, the parameter /3 must take the value 

4 


/3 = 




,3/2 


8 M 


(15) 


Thus the metric function is 


/(r) = 1 + - 
nr 


, 32M^r^ + SnMq^r + 

In--o-vr “ 2 tan 

32M^r^ - SnMq^r + n^q^ 


-T 


8Mr 


V 


+ 11+2 tan ^ ( 1 — 


8Mr\ 

7rg2 y 


(16) 


Expanding the metric function near r = 0, one can find that there is indeed a de Sitter core, 

= (17) 

3 [7r^q°) 

And expanding the metric in the r —)■ oo limit, we can obtain 

2M 

/(r) = l-+ % + 0(r)6. (18) 


Obviously, it recovers the RN black hole in the r —)• oo limit. 

There are three variables,(r, q, M) in the metric function. We can introduce a new coordinate 
X = r/M and a new charge Q = M/q. The parameters M, q can be cancelled in the metric function, 
we obtain 


/(x) = 1 + 


TTX 


32(3^x^ + 87r(5^x + TT^ _i 

In-- K -+ 2 tan 

32Q^x^ - 8nQ^x + 


1 - 


8Q‘^x 


TT 


— 2 tan 


-1 


8Q‘^x 


+ 1 


TT 


• (19) 


The horizons correspond to the locations where /(x) = 0. The equation is too complicated to 
express the roots analytically. As is shown in Fig.l, there can be one or two horizons depending 
on the value of Q. Particularly, for certain choices of Q, there can be no horizon, leading to a 
particle-like solution. For fixed x, if we decrease the value of Q, these two horizons of the regular 
black hole come closer and they meet at the critical value of Q = Qe- In this time, the black 
hole is extremal. Numerically, one has Xg = 0.9526 and Qe = 0.9893. This means that only when 
M > 0.9893(/ the horizons may exist. 


IV. WEAK ENERGY CONDITION 

The weak energy condition(WEC) corresponds to the statement T^yU^u'^ > 0, which means that 
the energy density of any matter distribution as measured by any local observer in spacetime must 
be nonnegative and dominates over the pressure[26]. Except for some special cases or matter fields, 
the WEC is fulfilled by general matter fields. Therefore we expect that a physically acceptable 
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f(x) 



FIG. 1: The metric function f{x) as function of x. The dashed line, solid(blue, thick) line, and solid(red) 
line correspond to Q = 0.8, 0.9893, 1.2, respectively. 

NED model cannot violate the WEC. According to the metric ansatz and Eq.(ll), WEC can be 
expressed equivalently as the following inequalities [27]: 

1 dm{r) 

r‘^ dr ~ ' 

From Eq.(14), it can be easily derived that 

dm{r) 

dr ^j3q^ + 2r^ ’ 

and 

2 dm{r) d^m{r) 2q^r^ 

r dr dr"^ ( 4 / 3^2 ^4)2' 

Clearly, the NED model we considered satisfies the WEC. 


2 dm{r) 
r dr 


> 


d^m{r) 

dr^ 


q2r2 


( 20 ) 




V. THERMODYNAMICS 


We now would like to study the thermodynamic properties of the regular black hole we obtained. 
The Hawking temperature of the black hole can be calculated using the formula 


Th = — = - 

27r dvr y/—gugrr 


r=r+ 




(23) 


r=rj^ 


where k is the surface gravity. In fact, it is not necessary to substitute the expression m(r) in 
Eq.(14) into the above equation. According to Eq.(ll), we know that m(r_|_) = r_|_/2 at the event 
horizon. Moreover, using the field equation, Eq.(12), one can easily obtain 

1 r_|_£ 


Th = 


47rr_i 


Svr 


(24) 
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We also use the new variables x, Q to express the Hawking temperature, which is 

^ 1024Q8xi - W24Q^xl + vr^ 

^ 40967rgQ^x^ + 4:TT^qQx+ 

In Fig.2, we show the temperatures of the regular black hole. It is shown that the black hole will 
never evaporate completely and radiate until the temperature becomes zero at some critical value 
of the event horizon. After that, the temperature is negative, which is physically unacceptable^. 
The zero-temperature black hole is in fact the extremal one. Here one thing should be clarified. For 
the metric function depicted in Fig.l the extremal case corresponds to fixed x = Xe and Q = Qe 
because the metric is the function of two variables x, Q only. While the Hawking temperature 
depicted in Fig.2 is dependent on three variables x+, q, Q. Thus, choosing different q, Q, the 
black hole always has the extremal case. 






FIG. 2: qTh as functions of x+ for (5 = 1, (5 = 1.5. The subfigure on the RHS is a magnification of the 
subfigure on the LHS in the range x G [0,1]. 


Now we should discuss other thermodynamic quantities of the regular black hole, such as internal 
energy and entropy and heat capacity. Entropy, Komar energy, Smarr formula, phase transition, 
etc. of regular black holes have been studied extensively [28-38]. However, the definitions of the 
internal energy and the entropy of regular black holes has some controversies [39]. Because our 
black hole solution is derived in general relativity, we insist that its entropy should satisfy the 
Bekenstein-Hawking area law, namely S = A^/A = vrr^. Generally, for regular black holes, if one 
takes the black hole entropy as the area-law form, the black hole mass M is no longer the internal 
energy of the system and thus the usual first law of black hole thermodynamics may violate [39]. 
However, black hole is a thermodynamic system, the first law of thermodynamics must be satisfied. 
Therefore, there must be some internal energy U which fulfill the first law: dU = TdS + (j)dq, where 

^ Although there are some physically acceptable systems which have negative temperature, for black holes the same 
explanation does not work. 
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q is the magnetic charge and cf) is the conjugate potential. To understand the local stability of the 
regular black hole, we can calculate its heat capacity. Although we do not know the detail form of 
the internal energy, we can still define the heat capacity at constant charge as 

_dU _j,dS _ 2-Kq^Q‘^xl {l024Q^xl + TT^) (l024Q^x^ - 1024^6^;^ + vr^) 

q - “ 1048576gi6a;8 _ 3145728^143,^ + 20487r4g83.4 ^ I0247r4g6a,^ + ^^8 ' 

(26) 



As is shown clearly in Fig.3 that Cq will approach zero as the regular black hole tends to the 
extremal one and diverge at the point where the temperature of the black hole takes the maximum. 
According to Davies’ viewpoint, this means that there will be a phase transition at Xj^ = Xm = 1-71 
with fixed <5 = 1 for the black hole. The Cq curve on the left of the point x = Xc, thus certainly the 
divergent point at smaller is meaningless because the temperature there is negative. The heat 
capacity is negative for larger black hole with > Xm and positive for smaller black hole with 
Xc < x+ < Xm- It means that the smaller black hole is thermodynamically stable and there may 
be a larger black hole/smaller black hole phase transition. At the critical point x^., a second-order 
phase transition happens. 

In addition, we can consider the thermodynamic quantities of the regular black hole from another 
perspective, namely the horizon thermodynamics proposed by Padmanabhan[40]. For a static, 
spherically symmetric spacetime the Einstein’s equations can be interpreted as the thermodynamic 
identity 

dE = TdS - PdV, (27) 

where T = Th, S = Ahj4:, E = y^A/i/lOvr, V = 47rr))/3 , and for our action P = —T^^(r/i)/87r. 

One cannot define heat capacity as C = dE/dT directly as is done in[41]. Firstly, due to 
the presence of matter fields, P ^ 0, thus the identity Eq.(27) cannot return to the usual first 
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law of black hole thermodynamics dM = TdS. Secondly, even in vacuum case, for many regular 
black holes dM ^ TdS with the area-law entropy [39]. Thirdly, in the framework of horizon 
thermodynamics one cannot define the heat capacity at constant volume Cy because that means 
constant r_|_. Therefore, we introduce the enthalpy H = E + PV of the system to define heat 
capacity at constant pressure as 



dH 


dH 

drp 

p 

dT 

P~ dT 



drp 

p 


According to Eq.(24), one can easily obtain 


Cp 


27rr^ (SvrPr^ -|- l) 
SnPr\ - 1 


(28) 


(29) 


Obviously, this result do not rely on the concrete black hole solution. It only depends on the 
gravitational theory under consideration. For general relativity, the heat capacity at constant 
pressure is given by Eq.(4), while for other theories of gravity it has different forms [42]. For our 
NED model. 


C _ -l2SM^q^ 

IGvr 10247rM^r;]L + ~ ’ 


(30) 



FIG. 4: The heat capacity at constant pressure Cp as function of r+ for different pressure. 

Obviously the behaviors of Cp and Cq are very different. This is because the two heat capacities 
are derived in two different frameworks of black hole thermodynamics. Cp does not diverge. 
However, there is a point after that Cp is positive which means that the black holes in this case is 
local stable. As we have noted above, this conclusion applies to all kinds of black holes derived in 
general relativity, not only our regular black hole. 
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VI. CONCLUDING REMARKS 

In this paper, we considered a model of nonlinear electrodynamics. We expect that the model is 
physically allowable because of three reasons. Firstly, the NED model returns back to the Maxwell 
electrodynamics in the weal field limit. Secondly, the NED model can produce the phenomenon 
of vacuum birefringence which is a kind of nonlinear effect originated from the nonlinearity of 
the NED. Thirdly, the electromagnetic field in the NED model satisfy the weak energy condition 
which means the energy density of the matter fields measured by a local observer is nonnegative. 
By coupling the NED to general relativity, we obtain a magnetically charged regular black hole only 

/ 3/2 \ ^ 

when the parameter /3 in the NED model satisfies the condition /3 = I 1 . We also shown that 
only when M > 0.9893g the regular black hole exists with one or two horizons. When M < 0.9893g 
we only have a particle-like solution. 

We also discuss the thermodynamic quantities and thermodynamic properties of the regular 
black hole. Erom Eig.2 one can see that the Hawking temperature becomes zero as the black 
hole tends to the extremal one. After that, the temperature is negative, which we think it is 
meaningless for black hole system. The regular black hole may have the similar problems to those 
we have studied in [39]. To discuss the local thermodynamic stability of the regular black hole, 
we calculate the heat capacity at constant magnetic charge Cq according to the corrected first law 
of black hole thermodynamics. It is shown that Cq is negative for larger black hole and positive 
for smaller one, which means that the smaller black hole is more stable. At the point where the 
temperature takes maximum, Cq will diverge. According to Davies’s viewpoint, a second-order 
phase transition happens there. 

In addition, we investigate the thermodynamics of the regular black hole according to Padman- 
abhan’s horizon thermodynamics. Dymnikova has ever employed this method to discuss thermody¬ 
namics of regular black holes[4I]. However, the heat capacity she defined is problematic. In fact, in 
the framework of horizon thermodynamics one can only define heat capacity at constant pressure 
Cp. Furthermore, the Cp is only dependent on the gravitational theories, not the concrete black 
hole solutions derived in the gravitational theories. In our model, the pressure is always negative 
so long as the mass or magnetic charge are nonzero. As is shown in Fig.4, the behavior of Cp is 
very different from that of Cq as one expected. 
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